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CHAPTER I. INTRODUCTION 


A. The Video Inertial Pointing System 

The video inertial pointing (VIP) system had its origin in the 
infrared (IR) astronomy program being conducted at the Ames Research 
Center, NASA, Moffett Field, California. In the VIP concept, a video 
sensor is used to provide data for the generation of pointing error 
signals and for a cathode ray tube (CRT) display of the star field. 

The CRT display significantly Improves the operator/astronomer inter- 
action with the remotely controlled telescope. Because balloon-borne, 
airborne, ‘or space-bome telescopes are typically gyro stabilized, the 
pointing error signals are used to update the basic gyro stabilization; 
hence the term video Inertial pointing. 

The system concept was originated by Mr. Charles D. Swift of 
Ames Research Center for use with a balloon-borne IR telescope. This 
early version of the VIP system [1] consisted of a silicon intensified 
target (SIT) vidicon boresighted to the telescope. On-board electronics 
were used to select the 10 brightest stars in the field of view of the 
SIT vidicon for transmission to the ground station and CRT display. In 
the ground station, the operator was able to select one guide star for 
automatic pointing. A joystick was used to surround the selected guide 
star with cross hairs that had an adjustable deadband. Telemetry was 
used to adjust the center position and deadband of the cross hairs in 
the on-board electronics. A coincidence detector determined if the 
star was in the cross hair deadband. If the star was within the dead- 
band, no command was sent to the gyro stabilization; If the star was 
outside the deadband, a fixed torquing current was sent to the appropri- 
ate gyro to recenter the star. During automatic pointing, the joystick 


1 



could be used to move the cross hair position and thereby introduce slew- 
ing commands to the telescope. Consequently, this system provided two-axis 
control (single guide star) and used a simple method of updating the 
gyro stabilization. Figure 1 is a block diagram of the airborne portion 
of this system, and figure 2 is a diagram of the ground station controls. 

Subsequent to the successful demonstration of this system on 
several balloon flights, the development of a more advanced version of 
the VIP system was started. The advanced system will include a charge 
coupled device (CCD) video sensor, three-axis control, and more sophis- 
ticated means of combining the pointing error signals with the gyro 
stabilization. 

B. Design Considerations 

The purpose of this work is to develop and evaluate a design method 
for the use of the video sensor and gyro data in updating the telescope’s 
gyro stabilization. The task is broken into two parts: (1) the use of 

the video sensor multi-star position measurements to determine the 
three-axis pointing errors; ‘and (2) the use of the pointing error signals 
combined with gyro data to update the gyro stabilization. The develop- 
ment of the multi-star processing and gyro filtering must take into 
account different control system analytical approaches and the perform- 
ance goals and hardware limitations of the VIP system. Performance con- 
siderations for a telescope pointing system include steady pointing state 
accuracy, pointing jitter or noise, and dynamic response to disturbance 
torques and slewing commands. 

The multi-star processing will directly affect the pointing accuracy 
and noise; the gyro filter will directly affect the pointing jitter, and 
will also affect the d 3 mamic response and the basic gyro stabilization. 

The gyro filter will not, however, directly affect the steady-state 
pointing accuracy, 2 
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Figure 2.- First generation VIP system - ground station. 
















The primary hardware limitation imposed by the VIP system is the 
use of a small digital computer to perform the multi-star processing 
and gyro filtering. Multi-star processing equations must be developed 
that do not require excessive computation time as this will adversely 
affect the overall system dynamics. Selection of an approach to the 
gyro filter must consider the computation time associated with the 
multi— star processing. Because the VIP system is intended to Interface 
with a telescope's basic gyro stabilization, an approach to the gyro 
filter design which does not require redesign of the gyro stabilization 
is desired. 

Although VIP may eventually be used with a variety of telescope 

gyro stabilization systems, the first use will be with a system that 

uses analog rate integrating gyros (RIGs) [2] . Consequently, filter 

designs developed in this work assume that use of the RIG-type of gyro 

and design parameters for this existing telescope gyro stabilization are 

used for the filter evaluation. 

■ Based on the limitations of the VIP microprocessor, the analytic 

approaches that can be used for the multi-star processing are somewhat 

limited. Consequently, small-angle, linearized equations are used and 

consideration is given to singularities and error performance. The 

resulting multi-star processing equations and selection criteria are 

believed to be unique. For the gyro filtering, a discrete filter, 

using the VIP microprocessor, is suggested due to the discrete nature 

of the video sensor star position measurements. Previous work [3] has 

considered the gyro -filter design independently of the design of the 

basic gyro stabilization and that is the approach used here. Additional 

justification for this approach is developed using a continuous’ filter. 

The use of the RIG integration in a continuous filter has been 
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previously suggested [3] and this is extended to the discrete filter. 
Since the video sensor and multi-star processing introduces time delays 
into the attitude error measurements, these must be included in the 
filter model. Design charts and analyses for a discrete filter with a 
similar application have been developed [4] but no consideration has 
been given to the inclusion of time delays or filter dynamics. A 
digital computer program, DISC [5] , based on the eigenvector decomposi- 
tion solution of the steady-state, optional filter problem is used to 
3-Tial.yze the different filter models. Both noise performance and 
dynamics are included in the filter analysis. 

C . Outline 

Chapter II describes the advanced VIP system hardware, operations, 
the CCD video sensor, and the operation of the CCD sensor. 

The development of the multi-star processing equations used to 
obtain the three-axis pointing error signals from the video sensor's 
multiple star position data is contained in Chapter III. Several 
<iiff6i^6nt methods are discussed and the resulting equations are 
examined for sensitivity and singularities. Based on the singularities, 
guide star and equation selection criteria are developed. The resulting 
equations are evaluated for accuracy performance using characteristics 
of the CCD video sensor. Performance improvement versus processing 
time is discussed and a summary of the multi-star processing approach 
concludes the chapter. 

Chapter IV describes the development of the gyro filter. A 
continuous filter is used to gain insight into and to justify an approach 

to the filter design which does not affect the gyro stabilization. A 

* •* 

steady-state discrete Kalman filter is developed for the VIP computer. 

Consideration of the system operation, video sensor, and multi-star 
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processing is used to determine what time delays must be included in 

the filter model. Six filters are developed that include different 

time delays and drift; they are analyzed for noise and pole locations 

for different combinations of sensor noise. Conclusions are drawn with 

respect to - selection of- sensor characteristics and resulting •'system" 

performance. A digital simulation is developed which models the gyro 

stabilization and telescope d303amics and includes a model of the video 

sensor and discrete filters. The simulation is used to verify the 

approach taken to the filter design and provides data that can be used 

to select filter characteristics for the VIP system hardware. 

> 

Chapter V presents results and conclusions; a discussion of areas 
in which further work is needed is included. 
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CHAPTER II. DESCRIPTION OF THE ADVANCED VIP SYSTEM 


A, System Description, 

The advanced VIP system [6] consists of three primary subsystems: 
(1) a CCD video sensor whose optical axis is coaligned with the tele- 
scope optical axis; (2) an electronics package, which includes a micro- 
processor and interface electronics; and (3) a control console with CRT 
display from which the system is operated. Figure 3 shows the VIP 
system and its primary subsystems. 

The video sensor used by VIP is a key element of the system. Two 
Important features of the VIP are dependent on the video sensor: the 
ability to provide three-axis control with a single sensor, and a CRT 
display of the telescope field of view for use by the operator. The 
CCD video sensor is an all solid-state device which uses interpolation 
techniques in order to obtain improved resolution. A microprocessor 
(in addition to the VIP microprocessor) in' the CCD sensor electronics 
contains algorithms for the digital interpolation. 

i 

The VIP electronics package has two main sections: a national 

semiconductor IMP-16 microprocessor, and a group of circuits that inter- 
face the microprocessor to the video sensor, telescope stabilization 
system, and VIP control console. Figure 4 is a block diagram of the VIP 
electronics. 

The microprocessor plays a central role in the VIP system. In 
addition to converting video sensor data into three-axis error signals 
via the multi-star processing algorithm, it performs all of the logic 
sequences necessary to initialize and track, it drives the CRT display, 
and it interrogates VIP control console switches to perform operator- 
controlled functions. The VIP microprocessor also communicates directly- 
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Figure 4.- Block diagram of the VIP microprocessor and interface electronics 






























with the video sensor’s microprocessor in order to effect data transfer 
as well as control the Integration time and threshold level of the 
sensor. 

The operator controls the position of the telescope line of sight 

the joystick on the VIP control panel and may observe a variety of 
®tar— field presentations on the CRT display. Switches controlling VIP 
operation, auto/manual select, fast/slow slew rate, and star— field dis- 
play options, surround the joystick controller. Threshold level and 
Integration period for the star tracker are selected by mpana of two 
thumbwheel switches on the operator’s left. 

Operation of a VIP— controlled telescope begins with coarse 
ac<juisition of the target star field. This is accomplished by slewing 
the telescope to a pre-computed orientation and making an initial align- 
ment with the joystick controller in the manual mode (joystick commands 
gyros directly). The star field, as viewed by the video sensor, is 
displayed on the CRT and the operator orients the telescope by com- 
paring the CRT picture with a star chart. Alternatively, a computer- 
stored star map of the desired viewing area may be presented. When the 
target has been identified, the automatic mode is selected by pushing 
the "initialize" button located on the end of the joystick. This 
commands the VIP microprocessor to store the x-y locations of all of 
the stars being tracked by the video sensor. The pointing .error in all 
three axes is calculated by comparing these current star positions from 
the video sensor with those stored in memory upon the "initialize" 

command. The calculated attitude error and the gyroscopically measured 

1 

attitude error are then combined in the microprocessor gyro filter to 
provide the updates to each of the three gyros. 
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The operator may continue to reposition the telescope optical axis 
during automatic operation using the joystick to change the stored star 
positions. A display option may be selected in which only the star 
field in its desired location is presented. Joystick activation moves 
this pattern on the CRT screen to whatever new position is desired. The 
dynamics of telescope motion are not seen (since only commanded star 
field position is displayed) thereby enabling very precise positioning. 
Simultaneous selection of current telescope position for the display 
indicates to the operator how close the telescope is to pointing in 
the desired direction. The CRT will show two star fields in this mode, 
a bright one for actual position and a dimmer field for the desired 
position. When the two fields are coincident, the telescope is on 
target. 

A digital indicator on the control console informs the operator of 
the number of stars currently ■ being used for tracking. Another display 
option permits viewing only these stars. As new stars enter or depart 
the field of view, VIP can display them without disturbing the tracking 
function, and it can keep the operator continually informed of the number 
and location of tracked stars. If, because of slewing or tracker param- 
eter changes, the ntimber or position of tracked stars is unsatisfactory, 
the operator may reinitialize the system by pressing the "Initialize" 
button again. 

B. CCD Sensor Description . 

The advanced VIP system is designed to use a sensor employing a 
charge coupled device (CCD) detector [7]. The face of the CCD detector 
is divided into a rectilinear pattern of discrete photosensitive ele- 
ments, each connected to readout registers also located on the front 
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surface. As photons strike the detector, the charge build up in each 
element is transferred to the readout registers after a selectable 
period of integration, the charge being directly proportional to the 
number of photons striking an element during the integration period. 

Since the location of each element is uniquely specified, it is possible 
to obtain the precise coordinates, within the resolution of the 
photosensitive-element matrix, of any stars being imaged on the surface 
of the CCD. 

The current state of the art for CCDs is a matrix of about 400 
elements square. Commercially available CCDs are about 200 elements 
on a side. Thus, even the most advanced units cannot obtain resolution 
better than 0.5 arc min for 3° FOV optics. However, it is possible 
through a process of interpolation to improve the resolution by an order 
of magnitude. 

The interpolation scheme requires that the output of a four by 
four matrix of elements be sampled and averaged each time the CCD is 
scanned. This array is centered on that element which first exceeds 
a pre-selected threshold during the initial scan of the device. Thus, 
for each star image, there is a corresponding -four by four matrix for 
which the charge "center of mass" is computed. It is possible, therefore, . 
to obtain computed star locations that are not centered on an element. 

The effect of this procedure is to Increase the resolution by artificially 
increasing the element density. 

The interpolation algorithm and the data formatting for use by the 
control system are complex, and their implementation requires extensive 
data handling. This is accomplished throjugh the_jus,e of a microprocessor 
built into the electronics of the CCD video sensor. Raw data from the 
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CCD are read at the end of each integration period by the microprocessor; 
the microprocessor then computes the locations of the 10 brightest stars 
in the field of view. Each star is assigned an identification number 
and the four pieces of data for each star (x, y, brightness, and i.d.) 
are transferred to the VIP computer when the star tracker microprocessor 
completes the data for all of the stars in a frame. Logic within the 
microprocessor automatically takes care of stars entering or leaving 
the field of view so that while it is always possible to display all 
stars in the field of view, the VIP controller only accepts data for 
valid stars. A valid star is one that has met several criteria, both in 
the video sensor and within the VIP software. Attitude information is 
generated using only valid stars. Initializations of the star tracker 
may be accomplished at any time at the option of the operator with a 
pushbutton located on the VIP joystick. 

The. integration time and threshold ar^ also controlled by the video 
sensor's microprocessor. Seven levels of both functions are selectable 
from thumbwheel switches on the VIP control panel. Figure 5 is a cutaway 
view of the star tracker that shows the major components, A modular 
approach has been taken to the optical and mechanical design so that the 
tracker may be easily modified to fit a variety of missions requirements. 
A Peltier-effect cooler is Included for CCD operation on the ground; 
operation of the cooler during high altitude or orbital flight is not 
required. 
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CHAPTER III. MULTI-STAR PROCESSING 


A. Introduction 

In order to determine the pointing error signals needed to point 
a telescope in three axes, measurements of the position of at least 
two known guide stars are required. A single guide star would provide 
pointing information for two axes perpendicular to the telescope’s line 
of sight, and gyro information, subject to drift, would have to be 
relied on to orient the telescope about its line of sight. • The video 
sensor used by the VIP system provides position measurements of several 
stars and these can be used to generate the complete three— axis point- 
ing error. Since only two stars are required to determine the three- 
axis pointing errors, other stars can be used to generate additional 
estimates of the pointing error, and averaging or statistical techniques 
can be used to reduce the effects of errors in the video sensor. 

There are several approaches that can be used to calculate the 
pointing error signals when the positions of two guide stars are mea- 
sured by the video sensor (video sensor boresighted to telescope line 
of sight) . One approach would use quaternions to relate the measured 
and desired pointing location; this method would be valid for large 
pointing errors, and for large angles between the guide stars and the 
video sensor. Alternatively, direction cosines could be used with the 
same result. Neither of these methods is desirable because the resulting 
equations involve trigonometric functions which are very time-consuming 
for a small digital computer such as the VIP microprocessor. Moreover, 
because the guide stars used by VIP are within a small angular distance 
of the telescope and the desired pointing location (video sensor has 
relatively narrow field of view), the large angle advantage provided by 
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these methods is not required. Consequently, small-angle approximations 
to coordinate transformations are used to develop the basic equations 
that provide the three-axis pointing error signals from position 
measurements of a guide star pair. 

B. Basic Equations 


If the position of a guide star is known with respect to the desired 
•telescope pointing position and if the actual telescope pointing position 
is in error by a small amount. 


i-^jj 


( 1 ) 


x! 1 0r -03^^ 

= -Or 1 0p 

-j-r L 6y -Op 1 _ 

where [ is defined to be the star direction in a vector basis aligned 
to the desired telescope pointing position and [ is the actual tele- 
scope- pointing position in the same vector basis, zj is aligned to the 
telescope line of sight and x^ and y^ are the transverse directions 
(yaw and pitch, respectively) . Since the video sensor measures only 
xj and y^, then 


= xj + Oryj - eyzj , y^ = -0rxj + y^. + Opz^ (2) 

If Xj and y^. are the desired location of the star, then the measured 
quantities (x,j., y^) and the desired quantities (xj, yj) can be subtracted' 

Ax 4 xj - xj , -Ay £ yj - y^ (3) 

to relate known quantities to the pointing errors (6r, 0 p, 6y) , Since 
the star direction can be considered a unit vector, zj = 1 and 

Ax = 8ryj - 0 y , Ay = -erxj. -b 0 p (4) 

Examination of these equations makes it apparent that one guide star 
provides two measurements related to the three unknowns, 0r, 0p, and 9y, 
and as expected at least two guide stars must be used to uniquely 

determine the three-axis pointing errors. 
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C» Processing Options 


The measurements of guide star position provided by the video 
sensor will be corrupted by errors caused by noise and nonlinearities 
in the sensor and by the errors resulting from the interpolation com- 
putations performed in the video sensor’s microprocessor. If the com- 
bination of errors is assumed to result in a Guassian distribution of 


star position measurement errors, then there are several options avail- 
able for using equations (4) for determining the pointing errors from 
several star position measurements. These options are discussed 
below. 


Option 1, Process all star measurement data simultaneously using 
a least mean squares approach (LMS) assuming prior knowledge of the 
error statistics associated with the star position measurements. For 


example, for three stars. 


and 


where 


Axi 




VI 

Ayi 


0r 


V2 

Ax2 

Ay2 

= A 

0p 


V3 

Vlf 

Axg 



V5 

_Ay3 




V6 


0r 

A 

®P 


PA%"1 


Axi 

Ayj 

Ax2 

Aj2 

Ax3 

Aya 


P 1 = aV^a, R = E[wT] , P = E[(0-e)(0-0)T] 


(5) 


( 6 ) 


(7) 


and 
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Option 2. Process star pair measurements using an LMS approach and 
then average the resulting pointing error estimates to obtain a final 
pointing error estimate. In a similar fashion to equations (5), (6), 
and (7), for two stars: 

(9) 


( 10 ) 

Assuming, independent errors for the x and y values of the star posi- 
tion measurements, the resulting equations for 0r, 6p, and 6y for 
each star pair are 







Ay^+^ 



( 11 ) 
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. . ; z ! o + 2 | o + 23 o . ila 

' Rl R2 R3 R4 

'xio ^ X 2 o^^ 


R4 y \R2 R4/\Ri R 2/ 


( 12 ) 


(13) 




Option 3. Process star pair data algebraically using equations (4) 
and then average the pointing error estimates obtained for each star 
pair to obtain the final pointing error. Since for two stars there are 
four equations in three unknowns, there are several independent ways to 
solve for 6 r, 0p, and 0y. The resulting equations are: 
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(15a) 


er 


Axj - Ax2 

^10 “ y2o 


or 


Ayi - Ay^ - 

^20 “ ^10 


(15b) 


ep 


xipAyg - x^pAy^ 
^10 ~ X20 


(16a) 


yipAyi + xiqAx 2 + y2pAyi - xipAxi 

y2o - yio 


(16b) 


X2qAX]^ y2pAy2 X2 qAx 2 + yipAy2 

yio - yzo 


(16c) 


_ yio^^2 - y2o^^i 

y2o - ^10 


(17a) 


yi0^y2 " + XgpAxi - yipAy^ 

XI p - X2P 


(17b) 


X2qAx 2 + y2pAyi + XiqAx 2 - y^pAy^ 
X20 - xio 


(17c) 


D. Comparison of Options 

The important items to be considered in comparing the three options 
are the number of arithmetic operations required and the pointing error 
performance. Option 1 can be eliminated at the ouset due to the large 
number of operations for the matrix multiplications and inversions when 
more than two guide stars are available. A comparison of the number 
of operations required for Options 2 and 3 can be made by considering 
the 0p equation (12) and the 0p equation (16b) or (16c). Table 1 
summarizes the number of algebraic operations for each equation, and it 
is easy to see that Option 3 offers considerable savings in computation 


21 



TABLE 1. COMPARISON OF OPERATIONS REQUIRED 
FOR STAR PAIR SOLUTION 


Operation 

0p 

Option 2 LMS 

Option 3 Algebraic 

Multiplication 

24 

4 

Division 

43 

1 

Addition 

Subtraction 

25 

4 


time. The number of operations for 0p assumes that the Ri’s for 

each star are known ahead of time but combinations such as [(1/Ri) + (I/R3) ] 

must be computed in real time after star pairs are selected. 

In addition to the number of arithmetric operations, pointing error 
performance of the two options must be considered. A simple comparison 
can be made by considering two stars on the x axis of the video sensor 
equally spaced about the origin (xiq = -X 2 q)* With this positioning of 
the stars, equation (16a) can be used for 6p and the variances of the 
error in 0p and 0p are 


E[4p] 


Rif + R2 

4 




R2R4 
R2 + R4 


From these, the ratio 


(18) 


' E[E0p] (R + R )2 

E[e? ] " ^R2R4 

0 P 

can be used to compare the relative performance of the two options. If 
R4 = R2» then (E[Sgp]/E[egp]) = 1 indicating no performance advantage 
for Option 2 when the star position measurement errors have equal 
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variances. When R 4 = 2 R 2 , then (E[£gp]/E[egp]) =1.125 and there Is 
large improvement, considering the large amount of extra 
computation time required. 

Based on this comparison. Option 3 appears to be the best approach. 
Examination of the equations for Option 3 reveals that there are singu- 
larities associated with each equation depending on the relative loca- 
tions of stars associated with each guide star pair. Consequently, 
consideration must be given to these singularities in the selection of 


which equation to use and in the selection of guide star pairs. 

E. Selection of Equations and Guide Stars 

Determination of which equations to use for 6 r, 6 p, 6 y and how 

to select guide star pairs can be made by considering the variances of 

the equations (15)- (17). If the errors associated with each Ax and Ay 

1 i 

measurement are assumed to be identically distributed and independent - 
with 


then, as an example. 


ElAx|] = E[Ay|] = 


E[£qj.] 


2g^ 

(yio - y2o)^ 


for equation (15a). Table 2 summarizes the 0r, 6 p, and 0y pointing 
error equations and their variances. 

The variances of the 0r pointing error equations, (15a) and (15b), 
depend on l/(yiQ“y 2 o)^ ~ respectively. A nondimen- 

sionalized plot of l/(xio-y 2 o)^ is shown in figure ' 6 .- .Although not 
shown, the first and third quadrants, and second and -'fourth quadrants 
are identical and there is clearly a discontinuity at the origin, A 
similar figure would hold for l/(yio “ yao^ The key feature of the 
diagram is the singularity when xio = X 20 or yio = y 20 the larger 
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TABLE 2. SUMMARY OF POINTING ERROR EQUATIONS 
E[Axi^] = E[Ayi^] = 



Attitude error equation 

E [Attitude error^]^ 

Singularities 


Axi - Axo 
— 1 


2cj2 

710=720 

Best performance 

i 

0r 

yio - 720 


<7lo - 720>^ 


7io yjo are 

opposite sign 


15b 

X20 ~ Xio 

2<j2 

(X20 - Xio)^ 

’'lO “X 20 

Best performance 
xjo and X 20 are 
opposite sign 


xioAy2-X2oAyi 
k(xio-X2o) ^ 

a^(xfo+x|o) 

XlO “X20 

Amplification of 
error when xjq 
and X 20 same 
sign 


k2(xjo -X2o)^ 


0p 

-710^71 +XioAx2 + y2oi7i -XjqAxi 

•16b 

o"i:(yL-^yin)-»-2x?j 

7l0 =720 

'Amplification of 
error when y^j 

k(y 2 o ” 7io)’ 

k^(720 -7io)^ 


and y 2 fl same 
sign 


X2qAxi - Vzo^Yz - ^ZO^^Z + 7 1 n^7. 

16c 

<^^[(ylo + yio>+24ol 

7l0 = 720 

Amplification of 
error when yjg 
and Yzo same 
sign 


. ^<yio-y2o> 

k^(720 -yio)^ 

t 


7ioiX2-y2flAxi 


<^^(710-^720) 

7l0 =720 

Amplification of 
error when y^g 
and y 2 o same 
sign 


k{y2o“7io) 


><'^(720 -710)^ 


6y 

710^^2 " +’'20^’'! - 7in*7i 

16b 

o2[(xJo+x2o)+2y^J 

X10 = X20 

Amplification of 
error when x^q 
and X 20 same 
sign 

k(xio - X2o) 

k^(xio-X2o)^ 



-X 20 AX 2 + y2oAyi +xioAxp - y^nAy^ 

16c 

a2[(x?Q + x|n)+2y2n] 

^10 = X20 

Amplification of 
error when xig 
and X 20 same 
sign 

a* ... 

k(x2o “ xjo) 

k2(xio -X2 o)2 



Assumes Ax^ and Ay^ are identically distributed independent random variables. 
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XiO 


Figure B.- Plot of l/(xio “ ^20)^ 






values when xiq and X20» or yig and 720 are the same sign. Consequently, 
when the x coordinates of the guide star pair are the same sign, equa- 
tion (15a) should be used, and when the y coordinates are the same sign, 
equation (15b) should be used for 6r. This selection is summarized in 
table 3. The case where two stars are in the same quadrant is eliminated 
by this selection strategy and results in a slight reduction in the 
probability of obtaining useful guide stars in a given field of view. 

This is not a serious problem and can be taken into account when design- 
ing the optics of the video sensor by considering lens size and field of 
view versus sensitivity and star densities. 

In a similar fashion, selection of the appropriate equation for 
0p or 0y can be determined. Since the basic form of the equations 
is the same for 0p and 0y, the selection rules for 0p can be used 
for 0y by observing the proper symmetry. Considering 0p, the vari- 
ance of equation (16a) depends on 



(xio - X2o)^ 


and a nondimensionalized plot of this equation is shown in figure 7. 
Again, the first and third, and second and fourth quadrants are identical 
and there is a discontinuity at the origin. The key features are the 
singularity when xjq = X 20 and -the larger errors resulting when 
xio and X 20 are the same sign. The variances of equations (16b) and 
(16c) for 0p can be rewritten as 


(16b) 


+ yip) 2g^xfo 
^^20 “ ^ 10 ^^ (^20 - yio )^ 


( 20 ) 
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•" TABLE 3. STAR PAIR AND EQUATION SELECTION CRITERIA 













^20 

A 


ho 

00 


Xl0^ + X20^ 
(Xio-X2o)^ 


-5 


-4 


-3 


-2 


-1 


.5 

.51 

.53 

.59 

.72 - 

-^1 
~ 4 

1 

~ 3 

1 

~ 2 

1 

" ^ 1 

1.6 

3.2 

8.5 

41 

CC 

.51 

.5 

.51 

.55 

.68 - 

1.9 

5 

25 

OO 

41 

.53 

.51 

.5 

.52 

.63 - 

2.5 

13 

OO 

25 

8.5 

.59 

.55 

.52 

.5 

.55 - 

5 

OO 

13 

5 

3.2 

.72 

"'68 

.63 

.55 

.5 - 

OO 

5 

2.5 

1.9 

1.6 


Xio 



I 


J> 





(16c) 


( 21 ) 


+ y?n) , 2g2x|p 

(y2o - yio)^ <y2o - yio^^ 

so that these equations have the same singularities when y = y and 

20 ■' 10 ■ 

the same amplification of error when and Y 2 o are the same sign. 

Consequently, when the x coordinates of the guide star pair are the same 
sign, equation (16b) should be used if ^ j^2ol* Alternatively, 

equation (16c) should be used if [x 2 ol < |xio|. If the y coor- 
dinates of the guide stars are the same sign, then (16a) should be 
used. The case where both guide stars are in the same quadrant is 
eliminated in a way sxmilar to that of 0r. Table 3 also summarizes 
the 0r, 0p, and 0y guide star and the equation selection process. 

F. Roll Sensitivity 

Although the selection process described in the previous section 
can be used, an additional criterion should be considered for the 
selection of guide star pairs for 0r. Comparing the variance of 
equation (15b) (0r) with (16a) (0p) we find that | 

E[e0rl 2 ’ ’ ' 

( 22 ) 

E[£0p] (xio - X2o)^ 

and if the xjq and X 2 q values are measured in degrees then 


^ 2(57.3)^ 

E[£0p] C^IO - ^20)^ 


Consequently, the roll pointing errors are much larger than the pitch 

and yaw errors unless the guide stars are far from the origin. This 

♦ 

result is not unique to the equations used, but is a result of the geom- 
etry of using stars close to the desired pointing direction to determine 

r 

the pointing error for the axis (0r) aligned to the line of sight. This 
suggests an additional criterion for selection of roll guide star pairs; 
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can 


select roll guide star pairs where one of the stars is outside of a 
radius with a pre-selected diameter. For this criterion, figure 8 
be used to determine an upper bound on the 0r errors as a function 
of this radius. 

G. Performance and System Tradeoffs 

Ultimate pointing performance of the system is a direct function 
of the video sensor errors, both steady state and noise, the multi- 
star processing, the gyro noise and gyro filtering. Video sensor noise 
sources include random events in the photoelectric processes, electronics 
noise, and quantization levels in the sensor's microprocessor. The num- 
ber of bits used in the sensor's microprocessor and in the electronics 
can be designed so that the primary noise source is the photoelectric 
process. Video sensor steady-state errors result from nonlinearities 
in the sensor and from the accuracy with which the sensor performs the 
interpolation calculations. 

The multi-star processing will improve both the noise and steady- 
state errors associated with the video sensor star position measure- 
ments. The reduction in errors results from the averaging and thus 
depends on the number of guide star pairs used. As an example of this, 
we can consider the CCD sensor. The largest error for the CCD sensor 
is the accuracy of the interpolation algorithms, and well spacing and 
non-uniformities limit the CCD interpolation accuracy to 10% of the • 

CCD array well size [7]. If the interpolation error is assumed to be 
uniformly distributed then 

E[Ax 2 or Ayp = ^ A o2 (24) 

where a is 0.1 times the well dimension in arc sec. Referring to 

table 2, and assuming guide star pairs that allow use of the "a" 
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equations for pitch and yaw 

E[0p^ or 0y2] < k crp^ or ay^ (25) 

since 


J. y.2. 

^10 - ^20 
(^10 - X2o)^ 


< 1 , 


^10 ^20 


(y 


10 


- 72o) 


2 — 


< 1 


(26) 


from figure 7 . If four guide star pairs are assumed then 




2 -, 0“^ 


N = 4 


(27) 


Similarly, if guide star pairs with one star of each pair are assumed to 
be at least 2“ from the origin 




(28) 


from figure 6. If four guide star pairs are assumed then 


E[s 


0r 


AVE 


N-4 


(29) 


Table 4 summarizes the 1-cr roll, pitch, and yaw pointing errors result- 
ing from the interpolation error for a CCD array size of 400 x 4Q0 ele- 
ments. The relatively poor roll performance 'is a result of geometry, 
as explained earlier. 


TABLE 4. 1-CT POINTING ERRORS (ARC SEC) USING 
FOUR GUIDE STAR PAIRS 


Array size 
(field of view) 

a 


^f^AVE 

a 

p>y 

^P’^AVE 

400 X 400 
(2° x2“) 

1.8 

^44 

<22 

^.0 

40.5 
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The effect of the multi-star processing averaging is shown in 
table 4 and provides a 1/)M (N number of guide star pairs) improve- 
ment in the 1-a values of the pointing error signals. However, the 
time required to process the additional guide star pairs and to perform 
the averaging, increase proportional to N- Figure 9 shows these two 
effects, and it is clear that the computation time increases signifi- 
cantly faster than the errors improve. As discussed earlier, there 
are two compqnents of the overall system pointing error: steady-state 

and noise or jitter. The video sensor’s star position measurements 
contain both steady-state and noise errors and the averaging in the 
multi-star processing reduces both types of error equally. However j 
while the gyro filter does not affect steady-state pointing error of 
the system, one of its primary purposes is to combine the pointing 
^error signals from the multi-star processing with the gyro signal to 
obtain a low noise estimate of the pointing error signals. As dis- 
cussed in the gyro filtering section, the time required to perform the 
multi— star processing Increases the delay between actual star position 
measurement and the time at which the resulting pointing error signals 
are used to update the gyro filter. It is shown that larger time delays 
increase the noise of the gyro filter’s pointing error estimate. Con- 
sequently, in a final system design, a careful tradeoff must be made 
between the reduction in steady-state errors provided by using addi- 
tional guide star pairs in the multi-star processing and the effect on 
noise performance of the gyro filter caused by the increased time delay. 
H. Summary of Multi-Star Processing Technique 

The technique developed for the multi-star processing using the 
linearized small angle equations can be summarized as follows: 
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COMPUTATION 
TIME 
TIMEN^N y 
TIME N = 1 X 


ERROR IMPROVEMENT 

ERROR N = 1 
ERROR N = N 


Figure 9.- 


2 3 
NUMBER OF GUIDE 
STAR PAIRS 

Comparison of computation time and error improvement versus number of. guide star pairs 



1. Select valid guide star pairs from those stars provided by the 
video sensor eliminating as pairs those stars that are in the same 
quadrant.. Additional criteria for roll is the elimination of potential 
guide star pairs unless one of the stars is farther from the origin than 
a pre-selected angle, 

2. For each guide star pair for pitch and yaw, select the appropri- 

I 

ate equation (table 3) and solve for 9p and 6y. For the roll guide 
star pairs, again select the appropriate equation and solve for 0r for 
each guide star pair. 

3. ' Average the - 0r, 9p, and 6y pointing errors obtained for each 

i 

guide star pair to obtain a final estimate of the three axis pointing 
errors . 


35 



CHAPTER IV. GYRO FILTERING 


A. Introduction 

The multiple star measurements from the video sensor are used in 
< c 

the VIP microprocessor to derive the three-axis pointing 'error signals. 

These signals are then available for use in the overall telescope con- 

. r 

trol system. For a typical telescope pointing system, both gyro and 

star sensor measurements are available. A typical system will include 

a basic gyro stabilized gimbal using a high performance gyro with the 

c 

star sensor measurements used to update the gyro. Since the system 
that will be used for the VIP system demonstrations and development 
uses analog rate integrating gyros (RIGs) , the goal is to develop a 
gyro filter that will combine the three-axis pointing error signals 
available in the VIP microprocessor with the RIGs to provide the best 
possible performance. 

The approach to be taken is to develop a discrete steady-state 
Kalman filter that combines the attitude error signals with the data 
from the RIG in the VIP microprocessor to provide update signals to 

s 

the RIG and the stabilization system. Because it is intended that 
the VIP system be able to interface with different telescope pointing 
systems (i.e., the basic gyro stabilized gimbal), the filter design and 
implementation must not affect the basic gyro stabilization. The design 
of the filter must take into account the processing delays associated 
with the video sensors processing and Interpolation and the multi-star 
processing time in the VIP microprocessor. 

B. Design Considerations (Continuous Filter) 

Although the attitude measurements from the video sensor and the 
multi-star processing ar'e available at discrete intervals and although 
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a discrete filter will be developed, some of the key aspects of the 
design, approach can be .most easily illustrated by considering a con- 
tinuous filter. Figure 10 represents a basic gyro stabilized girabal. 
If continuous measurements of 9 are available, then a simple filter 
can be developed to combine the attitude measurements with the RIG to 
provide a best estimate of 9 for use in the gyro stabilization loop- 
If 


■ 9 = 0), o)g = 0) + (Ogjj 

where ojgjj = gyro rate noise, and Wg = gyro rate signal, then 

0 = 0)g - 0)gjj 

and if the attitude measurement 


(30) 


(31) 


y = 0 + On (32) 

where 0 ^ telescope attitude and k attitude measurement noise, 
then a simple filter can be developed as 

0 = o)g + Ki(z - 9) (33) 

where Kj can be selected based on the desired noise and dynamic 
response. A block diagram of this filter is shown in figure 11 and it 
is easily seen that the filter can be incorporated into the gyro stabili- 
zation loop as shown in figure 12, where the RIG provides the integration 
required in the filter. This result has been previously noted [3] . 

Some important properties of this system can be derived considering 
simple RIG and telescope dynamics as shown in figure 13, where Kj has 
been generalized to G(S) . The first item of interest is the disturbance 
response, 0/Tjj, and it is easily shown that 


A = 1 

T^ s2 + H(s) 


( 34 ) 
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Figure 12.- Gyro stabilizafiion with ^filter. 
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MEASUREMENT) | + 


©b (JOYSTICK BIAS) 

Figure 13.- System with generalized filter. 







showing that the disturbance response is unaffected by the filter 

dynamics. For response of the telescope to a change in bias pointing 
position. 


- G(S)H(s) 

0s [s + G(S)][s + H(s)] 


(35) 


indxcatrng dependence on both the gyro stabilization and filter dynamics 
as expected and showing a separation effect, where the system closed 
loop poles are the poles of the filter and the poles of the basic gyro 
stabilization obtained independently. 

Considering the noise performance of the system, 

^ G(S)H(s) 

0N [s + G(S)j[s + H(s)] 

and 


9 H(s) 

“gN [s + G(S)] [s.+ H(s)] 


(37) 


If the performance of the filter, removed from the system, is considered, 

9 ^ 1 

' ^gN ® ®(S) , (38) 

and 


J_ ^ G(S) 

9n s + G(S) (39) 

If the gain and bandwidth of H(s), the gyro stabilization compensation, 
are sufficiently high as is usually the case 


H(s) „ 
s2 + H(s) “ 


over the frequency range of interest. This implies that the overall 
noise performance of the system will be close to the noise level 
predicted by the filter analysis. 
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The conclusions drawn from this analysis are: 

1. The disturbance torque response of the system is determined 
by the gyro stabilization dynamics and is not affected by the filter 
dynamics , 

2. The system response to a bias change in attitude is determined 
by both the filter and gyro stabilization dynamics. 

3. The overall noise response of the system (gyro and video sensor 
noise) will be closely predicted by the noise given by the filter 
analysis. 

These conclusions suggest that the filter can be designed independently 
of the basic gyro stabilization, and this approach will be used in the 
following sections. Although a continuous system was used for the dis- 
cussion, the general results carry over to the discrete filter, as long 
as the basic structure used in the analysis applies. Furthermore, the 
conclusions are valid only insofar as the filter design adequately 
reflects the actual dynamics of the system. This will be an important 

j 

consideration when the time delays associated with the video sensor and 
multi-star processing are considered, and these time delays will have 
to adequately modeled in the filter. ■ 

C. Discrete Filter 

t 

Since the attitude error signals are available at discrete instants 
rather than continuously, and since a digital computer is available, a 
discrete filter is suggested. From equations (30), a discrete model 
of the system is 


= 6(ti) +. (tOg - t0gjj)dt 

y(ti+i) = e(ti+i) + 0N(ti+i) 


(40) 
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and a discrete "Kalman" filter for this system is 


= ecti)*^ + ojg dt 

0(ti+l)+ = 9(ti+i)" + Ki[y(ti+i) - §(ti+i)“] 


( 41 ) 


Although the attitude error measurements are available only at discrete 

times, the gyro operates in continuous fashion, and it is desirable to 

implement the discrete filter making use of the gyro integration but 

not disturbing the basic gyro stabilization. Figure 14 is a block 

diagram of an implementation of the discrete filter that accomplishes 

this. The integration of the gyro is used to provide f^i+1 m dt 

^i ® 

and the output of the gyro plus the update is the best estimate of the 
attitude error. The gyro output plus the last update provides the best 
estimate of the attitude error in between updates. The discrete Kalman, 
filter can be used taking advantage of the gyro integration in a way 
similar to that of the continuous filter, and without breaking the basic 
high bandwidth gyro stabilization loop. Thus, analysis techniques avail- 
able for the discrete Kalman filter can be used to obtain noise and 
dynamic performance and this filter can be properly mechanized. 


More complex filters, which include gyro drift and time delays, 
can now be developed and analyzed. Time varying Kalman gains will not 
be considered because their computation would place an additional time 
burden on the VIP microprocessor, and constant gains should provide 
adequate response times. 

Gyro dynamics will not be included in the filter models, because 
typical RIGs used in gyro stabilization loops have fast response, and 
the gyro dynamics will have little effect on the filter performance. 
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Figure 1’4.- Implementation of discrete filter. 






D. Time Delays 

Before proceeding to the design and analysis of the digital filters, 
the time delays to be represented in the filters must be defined. These 
time delays result from the interpolation performed by the video sensor, 
and from the multi-star processing performed in the "VIP microprocessor. 

The video sensor has a selectable integration time ranging from 
approximately 100 msec to 1 sec. Larger integration times provide 
increased sensitivity and allow the use of dimmer stars- After an inte- 
gration cycle is complete, data are read from the CCD array and processed 
in the sensor's microprocessor to provide the interpolated position of 
the 10 brightest stars above a selected threshold. The Interpolation 
process takes place during the next integration cycle; this operation 
is illustrated in figure 15. The minimum integration time for the video 

I 

sensor is set by the time required for the data transfers and interpolation. 

The VIP microprocessor takes the multi-star position data from the 
video sensor; it then performs the multi— star processing and the gyro 
filter equations, outputs data to a CRT display and gyros, and interro- 
gates the control panel for commands. The VIP microprocessor operates 
asynchronously with the video sensor. When multi-star position measure- 
ments are available, the VIP microprocessor stops other activities and 
enters the data into memory. The most recent star position data are 
used in the multi-star processing. Upon completion of the multi-star 
processing, the gyro update equations are completed. 

Because outputting the CRT display data and interrogating the control 
panel are interspersed with the multi-star processing, the multi-star 
processing requires almost the full cycle time of the VIP microprocessor. 
This cycle time depends on both the number of stars processed and on the 
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number of display options, and is generally constant once a particular 
display mode and number of guide stars have been selected. The cycle 
time of the VIP microprocessor will range between 300 and 500 msec. It 
is assumed that the time required for the gyro filter update is a very 
small portion of this time. The asynchronous operation of the two 
microprocessors is illustrated in figure 16. The result of this as3m- 
chronous operation is the occasional loss of star position data when 
the video sensor integration time is less than the VIP microprocessor 
time, and the occasional use of old data when the integratidn time is 
longer than the VIP microprocessor cycle time. 

The ultimate time delay associated with the pointing error signals 
used to update the gyro filter in the VIP microprocessor depends on 
several things. One delay that is constant from cycle to cycle (VIP 
microprocessor) is the multi-star processing time which can be modeled 
as a unit delay of the VIP microprocessor .cycle time. The other 
principal delay depends on when the star position data have been read 
into the -VIP microprocessor,. This delay- is variable from cycle to 
cycle due to the asynchronous operation of the two microprocessors and 
will vary from the video sensor's interpolation time up to the full 
video sensor integration time; this is shown in figure 16. This does 
not count' the time during which the star data to be used are being 
Integrated in the video sensor; and this integration process is not 
modeled here. 

The actual time delay -will be variable from cycle to cycle (VIP 
microprocessor) and will range approximately from Tyjp to Ty^jp + -Tysi 
where TyiP c the VIP microprocessor cycle time and Tygj = the selected 
video sensor integration time. The approach used will be to model the 
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Figure 16.- Asynchronous operation 


time delay as unit delays of the VIP microprocessor, because the time 
delay is variable from cycle to cycle. For the case where Tysi < TviP> 
a filter with one-unit delay should be satisfactory, and when 
T-^si > Tvip, two-unit delays will probably be required. A more complex 
approach would be to use electronic circuits separate from the VIP 
microprocessor to determine when the video sensor data were available, 
and thus determine the actual time delay for each cycle. These data 
could be used at each gyro filter update to select a filter model or 
gains for the exact delay. It is not thought that this complexity is 
needed, and the simulation results given later in this chapter verify 
the approach taken. 

E. Filter Models 

In the approach used here for the filter design and analysis, we 
start with the simplest filter and examine successively more complex 
filters. This provides better insight into the changes in filter noise 
performance and dynamics as a function of the complexity, and will 
ultimately allow a better understanding of how to select filters for 
the actual hardware. Gyro drift will be included in the filter models 
along with zero unit delays, or one- or two-unit delays. The integra- 
tion process of the video sensor and the gyro high frequency dynamics 
are not included. 

With these assumptions, there are six filter models of interest; 
they are listed below with the simplest model listed first. 

1. One-State Filter 

0 Attitude 

2. Two-State Filter 

0 Attitude 

D Gyro Drift 
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3. 


Two-State Filter 


0 Attitude 

0g Delayed Attitude (Unit Delay) 

4. Three-State Filter 

9 Attitude 

D Gyro Drift 

9g Delayed Attitude (Unit Delay) 

5. Three-State Filter 

9 Attitude 

9gj Delayed Attitude (Unit Delay) 

0g2 I5elayed Attitude (Unit Delay) 

6. Four-State Filter 

9 Attitude 

D Gyro Drift 

6gj^ Delayed Attitude (Unit Delay) 

0g2 Delayed Attitude (Unit Delay) 

For each filter, DISC [S-] will be used to analyze filter noise and 
dynamic performance, and, therefore, the state models must be developed. 
For use in the actual VIP hardware, the implementation of each filter 
similar to figure 14 must also be developed. The input to DISC is in 
the form 


“ ^9(ti) + riu(ti) + r2u(ti) 

y(ti) = H0(ti) + v(t^) 


(42) 


For filters that do not include drift, o)(ti) ->• N(0,Qj) where depends 

on the time between updates and is the integral of the gyro rate noise. 

If Ugjj N(0,Qg), then Q^j = TQg where T is the time between updates. 
For filters that Include drift. 
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where Qg is. the gyro rate noise as before and Ditj is the drift rate 
noise [4] . 

The unit time delays are included as shown below for filter 3t 


Model: 


6(ti+i) ~ 0 (ti) + (ug-o)gjj)dt 


= 0(ti) 


y(ti+i) = 0 s(ti+i) +v(t£+j) 


Filter: 


ecti+p- = 

9(ti+i)+ = 
9s(ti+i)+ = 


6(ti)^ + tOg dt 

0(ti)+ 

9(ti+i)“+Ki [y(ti+i) -§sCti+i)"] 

e(ti+i)“+K 2 ^[y(ti+i) - 0s(ti+i)-] 


Table 5 shows the DISC models for the, filters and figures 17 
through 21 show' the hardware implementations. 

F, _ Observability 

The filters' which include the unit time delays have a singular 


matrix. Consider filter ,3 ( 9 , 80 ) 


^1 o'! 

-1 oj 


which is clearly singular due to the zero column. For processes defined 


X = Fx 


where F is finite, the resulting state transition matrix ^ is 


nonsingular and 
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TABLE 5. DISC MODELS 
= P 


Filter 


T2 

H 

Wti 


X 

1(6) 

[1] 

[1] 

[1] 

[Qd] 

[R] 

[6] 

2(0,D). 

r 1 -T 
Lo 1 

■1 


'1 0" 
.0 1_ 


[1 0] 


Qd+|Dj^T 3 -|dkT2" 
_-|dnT2 DkT _ 


[R] 


a 

3 ( 6 , 0 g) 

G S] 




[0 1] 

[Qd] 

[R] 


[:.] 

4 ( 6 ,D, 0 g) 


Pi -T O' 
0 10 
Ll 0 0. 



^1 0l 
0 1 
^0 oj 


[0 0 1] 

> 


Qd+foNT^ --|dnT2” 
__--|dnT2 DkT _ 


[R] 


— 1 

e 

D 

_®s_ 

5(6,032,632) 


fl 0 o' 
10 0 
Lo 1 0 _ 



3 

[0 0 1] 

■ [Qdl 

[R] 


■ 

6(e,D;0si>6s2) 

L_ 

"l -T 0 
0 10 
10 0 
0 10 

0 “ 

0 

0 

0 



"1 0 ~ 
0 1 
0 0 
_0 0^ 


[0 0 0 1] 


Qd+-|%T 3 

- DjjT 


[R] 

. 

■ 0 "1 

D 

0S1 

_0S2j 


















4> 



SAMPLE 
AND HOLD 
















Figure 18.- Implementation of filter 3 (0,6g). 
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Figure 20.- Implementation of filter 5 (e,0gj,0g2). 
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<i>(t2. ti) = (j) Hti, t2) (48) 

Therefore, x<t 2 ) can be found given x(ti), or x(ti) can be found given 
x(t 2 ). For our situation, the model is sufficient to find x(t 2 ), given 
x(tj) by delaying, but the model is not sufficient to find x(ti), given 
x(t 2 >; that is, (J>"^(t 2 , tj:) does not exist. Considering the error 
equation for the filter 

• • e(ti+i) = - KH<i>]e(ti) (49) 

where 


[<|> - KH(j.] = 

and the characteristic equation is 


1-Ki 0 

I-K 2 0 

Z(Z + Ki - 1). 


Clearly 


the error equation cannot be arbitrarily selected since K 2 
effect and one root is always at the origin. The system is 
however, since the rank [8] of 


• (50) 

the roots of 
has no 
observable. 


[hT I = 



I 


(51) 


is two. 

The singular ^ matrix does present a problem for the use of DISC, 
because DISC calculates <[) ^ [5] .■ This problem' can be eliminated by 
introducing a small' coupling term in iji. 




" 1 

.1 ej 


(52)' 


Results from DISC for-, filter 3, using different' values of e, are shoxm 
in table 6 and show the Insensitivity to e. As e gets smaller, one 
pole gets- closer to- the origin. For values of e below 0.00001, the 
resulting M and P matrices begin to get unsymmetrical . It is inter- 
esting that K 2 = Ki which is apparently important for minimization of 
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TABLE 6. DISC RESULTS (FILTER 3) VERSUS COUPLING TEEM (e); 

Qd = R = 1 


e 

Mg 

P0 

Kl 

K 2 

Pole 1 

Pole 2 

0.01 

2.6 

1.6 

0.62 

0.62 

0.37 

0.01 

0.001 

2.6 

1.6 

0.62 

0.62 

0.38 

0.001 

0.00001 

2.6 

1.6 

0.62 

0.62 

0.38 

0.00001 


the filter noise response, but K 2 has no apparent effect on the pole 
locations. Also K 2 is between 0 and 1 since if K 2 '= 0 

^ 0s(ti+i)‘*' = 9s(fci+i)“ = e(ti)+ - ■ (53) 

which ignores the current measurement, and if K 2 = 1 

SsCti+i)'*' = y(ti+i) (54) 

which ignores past measurements. 




e(t:i) 


(55) 


This approach can be further justified by considering the error 
equation for filter 3 with the coupling term, 

1 - Ki -sKi* 

1 “ K 2 -sK2 

and the characteristic equation is 

+ (Ki - 1 + eK2)Z - e(K2 - Ki) (56) 

As e gets small, the equation reduces to the characteristic equation 
obtained without the coupling term, and for small e the roots are 
approximately Zj = 1 - Ki and Z 2 = 0. 

The approach to be taken, then, will be to use DISC for the filter 
analysis with small coupling terms for those filters that include the 
•unit time delays. For filters 3 and 4, which include one delay, 

£ - 0.00001 is used, and for filters 5 and 6, which have two delays, an 
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e is required for each delay and s ~ 0.01 is used. A higher value 
of s is required when the filter includes two delays because 
terms result in the calculations and DISC has problems at correspondingly 
higher values of e. This approach allows the use of DISC for the optimal 
^^ 2 ,lysis which provides a filter optimized for minimum noise 
response. The filters that include the unit delays do not allow arbi- 
trary selection of filter poles. But the DISC analysis is not oriented 
to this design approach and the poles, which cannot be adjusted, are at 
or very close to Z = 0, and Z = 0 is associated with deadbeat response; 
that is, a response that decays in one cycle. 

G. Filter Analysis 

The purpose in this section is to analyze each of the six filters 
with DISC for various combinations of gyro and video sensor noise to 
obtain the filter noise performance (primarily 0) and dynamics. As 
previously stated, the system noise perfoi;mance will closely approxi- 
mate the filter noise performance (9), if other sources, such as random 
disturbances and torque motor stiction, are neglected. The filter 
dynamics will play a large role in determining the overall closed loop 
response to step changes in the command attitude or joystick bias. 

Each of the six filters is analyzed so that an intuitive feel for 
the effects of additional time delays and the inclusion of gyro drift 


can be gained. Appendix A contains a summary of the DISC runs made 
for each of the six filters. Table 7 shows the combinations of 
and R used for each filter. For filter 2, the values of 

Ud + J - I %t2] 


to(t^) 



%T 


(57) 
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TABLE 7. COMBINATIONS OF AND R FOR FILTER ANALYSIS 


Qd 

0.1 

1 

10 

100 


0.1 

0.1 

0.1 

0.1 

R 

1 

1 

1 

1 


10 

10 

10 

10 


100 

100 

100 

100 


Q(j = (arc sec/sec) ^ 
R k (arc sec)^ 


are tabulated in table 8 . Gyro drift rate noise, Djj, must be included 
so that all filter ^ates are driven by some disturbance noise. ^The 
DISC analysis for filters 4 and 6 was performed only with Djj = 1 
because the results from filter 2 with Djj - 0.1 indicated poor role 
locations and slow filter dynamics. 

Results for the one— state filter can be examined first since the 
more complex filters will relate to these results. It should be 
remembered that M 0 is the covariance of the attitude estimate just 
prior to a filter update, and that P 0 is the covariance just after a 
filter update. The total system (filter plus gyro stabiliation) noise 
response depends on both M 0 and P 0 since the gyro is the primary 
reference between filter updates, and M 0 is closely related to the 
gyro rate noise. 

For presentation of the data, M 0 and P 0 can be normalized with 
respect to R; .figure 22 is a plot of M 0 /R and P 0 /R for various 
Qji/R. For large Q^/R, Mq approaches Q^j, and P 0 approaches R 
as expected. For Q^/R small. Mg approaches Pg in between and 
R. Figure 23 is a plot of Mg/Pg versus Q<j/R, and illustrates these 
conclusions. Figure 24 is a plot of pole location versus Q<j/R; for 
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POLE 



large Qjj/R, the pole is close to Z ~ 0 indicating rapid decay of 
initial condition errors, and for Q^/R small, the pole is close to 
Z = 1 indicating slow response. 

From these data, it can be concluded that Q(j/R should be 
selected so that 


0.1 



10 


This ' is because of the following : 

(1) If Q{j/F. > 10, pole location will be good but M 9 /P 0 will be 
large indicating that the filter is not performing much smoothing. 

Pq/R ~ 1, that is, at the filter updates the filter follows the video 
sensor data. Since M 9 /P 0 is large, system noise performance will not 
have been improved a great deal by the filter. 

\ 

(2) If Q(j/^ 0.1, M 0 /P 0 will be close to 1, but pole location 

will be close to Z = 1 and dynamic perforaance will be poor. 

For -filter 2 ( 8 ,D), the results are similar to those for filter 1 
with the exception of the additional pole and somewhat larger M 0 and 
P 0 . The poles generally consist of two real poles, one of which depends 
on Q(i/S- and is near or equal to the pole of filter 1 , and the drift 
estimator pole which depends on as shown in figure 25. For 

large Qd/% the additional pole approaches Z = 1 and the filter 
response would be sluggish. For values of Qd/% small, M 0 and P 0 
are somewhat larger than for filter 1 and the effect is worse for 
larger values of T. For large, M 0 and P 0 are close to the 

values obtained from filter 1 . 

From filter 3 (6,0s) and filter 5 (e,egj, 0 g 2 )» which Include one- 
and two-unit delays, respectively, the primary effect is poorer noise 
performance; that is, larger values of M 9 and P 9 relative to filter 1. 
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For both filters, the gains are identical and equal to the gain 
obtained for filter 1. The response consists of one or two poles 
close to or at the origin plus the real pole obtained from filter 1 . 
Figure 26 is a comparison of M 0 for filters 3 and 5 versus filter 1, 
and figure 27 is a similar P 0 comparison. For small Q{j/F> the time 
delays do not have great effect while for large Q^/R, the effect is 
pronounced, particularly on Pg. 

For filter 4 (0,D,0g) and filter 6 (0,D,9g^ , 632 ) , which include 
both drift and one- and two-unit delays, respectively, the results are 
generally a composite of the simpler filters. Only the cases with 
drift noise Djj = 1 were used since it was found with filter 2 that 
Djj = 0.1 results in a pole close to Z = 1, even for values of 
Qd = 0.1 and 1.0. For filter 4, the pole locations are generally close 
to the pole locations for filter 2 , with the addition of a pole close 
to Z = 0 due to the unit delay. Noise performance (MgjPg) is gener- 
ally worse than that of filter 3 (0,0g) at low values of Qd/%» and 
approaches the values of filter 3 for large (see fig. 28) . A 

similar comparison can be made between filters 6 and 2 for pole loca- 
tion, and between filters 6 and 5 for noise performance (see fig. 29). 

It is important to note for the filters that contain gyro drift 
rate noise, Djj, that the large relative value of Djj used' for the DISC 
is not representative of the physical situation. For most medium to 
high performance rate integrating gyros, Djj will be much less than Q^. 
For filters 4 and 6 , this is representative of values of = 10 or 
100 where Djj = 1 , and for these cases, noise performance is relatively • 
unaffected by Dj^, and closely approximates the noise performance 
obtained for filters 3 and 5. The only problem mth large Qd/%, is 
the low drift gain and pole approaching Z = 1, which will cause poor 
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Figure 27.- P 


0 


comparison. 
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Figure 28.- Filter 4 (e,D,0g) comparison to filter 3 ( 0 , 63 ) 



Figure 29.- Filter 6 (6,D,9g-^,9g2) comparison to filter 5 (0,0g|^ ,0g2). 


dynamic response of the filter because the drift estimate pole appears 

r 

in the transfer function for 0. This problem can be eliminated, how- 
ever, by proper selection of the drift estimate gain. The comparison 
of filters 1 and 2 indicated that the drift estimate gain and associated 
pole location depend primarily on while the gain and pole loca- 

tion for the attitude estimate depends primarily on Q^/R- The implica- 
tion is that the drift estimate gain can be selected for a reasonable 
pole location independent of the selection of pole location of the atti- 
tude estimate. This will provide good dynamic response of the filter 
but there will be some increase in M 0 and P 9 for the attitude estimate; 
this can be seen by comparing the results in Appendix A for filters 1 
and 2 . 

The analysis of filter 1 suggested that 0.1 < Q^j/S. < 10 for good 
filter dynamic and noise performance. The analysis for the filters that 
include time delays showed that the ratio of P0 with time delays to Pg 
for filter 1 increased rapidly with Q^/R. Because Pg is relatively 
more important than Mg, a further reduction in the acceptable range of 
Q(j/R (fig. 30) where 0.1 < Q^/R <1, is suggested. 

H. Simulation and Results 

In order to verify the approach taken to the filter design and the 
results of the DISC analyses, a simulation of a typical gyro stabilized 
gimbal was developed for use on an IBM 360 digital computer. This simu- 
lation uses a simulation language, continuous system modeling program 
(CSMP) [9] , and models one of the gyro stabilized gimbals' of the AIROscope 
balloon gondola used by Ames Research Center [2] . Several of the filters 
analyzed by DISC were modeled, and various tests with sensor noise, 
torque disturbance, gyro drift and joystick commands were made to 
evaluate the system performance. 
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Figure 31 is a schematic block diagram of the simulation. The 
basic gyro stabilization is simulated in detail including; 

1. Gyro dynamics 

2. Gyro demodulator and prefilter 

3. Gyro stabilization compensation 

4. Torque motor with back EMF 

5. Saturation levels for: gyro; gyro demodulator and preamp; 

compensation amplifier; power amplifier; and torque motor 

The video sensor integration and interpolation are simulated by an 
integration and averaging process where 


and the video sensor output is 


(58) 




“•VSI 


(59) 


and ^VSI = video sensor integration time. 

The time delay associated with the interpolation is not simulated 
nor is the multi-star processing done in the VIP microprocessor. The 
multi-star processing time delay is simulated as one complete cycle of 
the VIP microprocessor. The asynchronous operation of the video sensor 
and VIP microprocessor is simulated. Appendix B contains a computer 
listing of the complete simulation with filter 4 (SsDjSg), and values 
for the simulation constants. 

From the discussion of the continuous filter, it was concluded that 
the filter could be designed separately from the gyro stabilization. If 
the discrete filter properly models the processing time delays, response 
to disturbance torques is determined by the gyro stabilization d 3 mamics. 
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JOYSTICK INPUT 


Simulation block diagram 























and system noise response (6) depends primarily on the filter noise 
performance. System response (0) to joystick commands depend on both 
the filter and gyro dynamics. 

Figure 32 shows the response of the basic gyro stabilization without 
video sensor feedback or the gyro filter to a step disturbance torque. 
Figures 33—35 show the response of the complete system to a step dis- 
turbance torque for filters with zero time delay and for those with one- 
and two— unit time delays, and with the video sensor integration time set 
at 0.25 sec and the VIP microprocessor cycle time set at 0.4 sec. The 
responses consist of the basic gyro stabilization dynamics in between 
the filter updates every 0.4 sec. The response with filter 3 (0,0g), 
which models one-unit delay, is rapid and most completely approximates 
the disturbance response of the gyro stabilization. This result was 
expected since one-unit delay most closely represents the actual delay. 

Figures 36-38 show the response to a step-disturbance torque with 
the video sensor integration time set at 0.9 sec and the VIP micro- 
processor cycle time set at 0.4 sec. The response with filter 5 
(0 , 0g2 > » which includes two-unit delays, is best, and the two-unit 
delays most closely represent the actual time delay. In this case, 

1'VSI “0-9 sec and Tyjp = 0.4 sec, the response with filter 1 (0) has a 
long term instability. Table 9 summarizes the disturbance torque 

TABLE 9. TIME REQUIRED TO SETTLE. WITHIN ±0.1 ARC SEC 
FOR Tq = 0.1 FT LB (FILTER GAINS = 0.62) 


System 

Tysi — 0.25 sec 
Tyip = 0.4 sec 

Tysi = 0.9 sec 
^VIP " 0.4 sec ' 

Filter 1 (0) 

4.4 sec 

Long term instability 

Filter 3 (0,0g) 

2.0 sec 

11.4 sec 

Filter 5 (0,0g2,0g2) 

16.5 sec 

2.3 sec 
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Figure 32.- Gyro stabillziation response. 
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Figure 35.- Filter 5 (e,0g^,0g2) 
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Figure 36.- Filter 1 (0) response 






sec 



Figure 37.- Filter 3 ( 6 , 63 ) response. 
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Figure 38.- Filter 5 (9,9g) T^j response. 






responses for the three filters and two video sensor integration times. 
The responses (filter 3 for Tysi =0.25 sec, and filter 5 for 
"^VSI =0.9 sec) are somewhat longer than the basic gyro stabilization 
response since we have chosen to model the variable time delays (asyn- 
chronous operation of the video sensor and VIP microprocessor) as unit 
delays in the filters. The responses are adequate, however, and verify 
the basic approach. 

For the joystick response, figures 39-41 show the system response 
to a 36 arc sec joystick command with Tysi = 0.25 sec and Tyip = 0.4 sec. 
As expected, filter 3 (0,0s) j which models one-unit delay, provides the 
best response. Figures 42-44 show similar responses for Tygj =0.9 sec 
and for this case, filter 5 ( 6 , 6 gj, 0 g 2 ) provides the best response. 

Table 10 summarizes the responses to the joystick commands. The results 
verify the choice of filter from the disturbance response for the two 
different video sensor Integration times. 

For the system noise response, filter 3 ( 6 , 63 ) was used with 
Tvsi = 0.25 sec and Tyjp = 0.4 sec. Attitude noise is added to the 
output of the simulated video sensor, as shown in the simulation block 
diagram, figure 31, and gyro rate noise is added to the input to the ' 
rate integrating gyro. To properly simulate the values of gyro rate 


TABLE 10. TIME REQUIRED TO SETTLE TO 36 ±1 ARC SEC 
JOYSTICK STEP AT T = 0 (FILTER GAINS = 0.62) 


Filter 

Tygj — 0.25 sec 
^VIP = 0.4 sec 

Tysi = 0.9 sec 
Tyip = 0.4 sec 

1 (0) 

6.3 sec 

Unstable 

3 .(0,0s) 

1.6 sec 

>20 sec 

5 (0>0sl>Ss2) 

18.2 sec 

5.8 sec 
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Figure 40.- Filter 3 (9,0g) joystick response. 
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Figure 42.- Filter 1 (0) joystick response. 
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Figure 44.- Filter 5 (6>6si>®s2^ joystick response. 


t 

noise (Q^) used in the DISC analysis, it was necessary to consider the 
integration time of the digital simulation [10] and x(k) = Q(j/At. 

Since At = 0.005 sec, the standard deviation of the gyro rate noise 
used in the simulation is 14.1 times the standard deviation of the gyro 
noise used in the DISC analysis. In table 11, the noise performance 
predicted by DISC and the actual simulation noise performance are 
compared. The standard deviation .of 0 was computed from the simulation 
results using a large number of data points at equally spaced time inter- 
vals and represents an approximation. Nevertheless the simulation 
results agree quite closely with the DISC analysis verifying the conclu- 
sion that the filter noise performance determines the overall system 
noise performance. 

In order to verify the basic approach taken to the filter design 
with the simulation, the drift term was not included in the filter 
models. To evaluate system performance w^th the drift estimator, filter 4 
(6,D,0g) is used with Tygj = 0.25 sec and Tyj-p = 0.4 sec. The system 
response to a drift step with drift estimator gains of -0.55 and -0.10 
is compared in table 12. The lower gain on the drift estimate results 
in slower filter dynamics in agreement with the DISC analysis, and the 
gain of -0.55 is used for the joystick command and disturbance torque 
responses to provide adequate dynamic response. Figure 45 shows the sys- 
tem response to a disturbance torque response and figure 46 shows the sys- 
tem response to a joystick command using filter 4 (0,D,0g) with the. drift 
estimator gain of -0.55. These responses compare directly to figures 34 
and 40, and in both cases the response time is significan,tly longer due 
to the effect of the drift estimator on the filter poles. 

Noise response with filter 4 (6,D,6g) is shown in table 13, and 
* / 

compares well with the values given in table 12 (no drift estimator) for 
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TABLE 11. NOISE EESPONSE: FILTER 3 ( 6 , 63 ); 

Tysi “ 0.25 SEC, Tyjp = 0.4 SEC 


Qd 

R 



K 

0 simulation 
standard 
deviation 

1 

1 

1.6 

1.3 

0.62 

1.3 

10 

1 

4.6 

3.3 

0.92 , 

3.1 

1 

10 

2.7 

1.9 

0.27 

1.8 

100 

1 

14.1 


0.99 

8.1 

1 

100 

3.4 

3.2 

0.1 

2.2 


TABLE 12. EESPONSE TO DRIFT STEP: FILTER 4 (0,D,0^); 

^ 0 

TVSI = 0*25 SEC, Tvip = 0.4 SEC; Ki = 0.62 


Drift gain 

Time to reach steady state ±10% 

-0.55 

6.8 sec 

-0.1 

12 . 8 sec 


TABLETS. NOISE RESPONSE: FILTER 4 (0,D,0g); 

Tygj = 0.25 SEC, Tyjp = 0.4 SEC 



R 

Ki 

K2 

0 simulation 
standard 
deviation 

1 

1 

0.62 

-0.55 

1.8 

10 

1 

0.92 

-0.55 

3.4 

1 

10 

0.27 

-0.55 

8.3 

1 

100 

0.1- 

-0.55 

60.5 
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Figure 45.- Filter 4 (0,D,0g) response. 
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Qd 1> K. 1, and Qd 10, R - 1. For the larger values of R/Qdj 
the noise response is much worse than the values in table 11. Since 
inclusion of the drift estimator in the filter causes longer responses 
to disturbance torques and to joystick commands, and poorer noise perform- 
ance, it may be desirable to use filters 'without the drift estimator 
during periods when low noise performance or rapid slewing with joystick 
commands is desired. 

I . Summary 

The approach taken to the gyro filtering design has been validated 
by the simulation. A discrete, steady-state Kalman filter is developed 
that uses the rate integrating gyro. Assumption of asynchronous opera- 
tion of the video sensor and VIP microprocessor lead to the use of unit 
time delays in the filter modeling the actual variable delay. DISC has 
been used to analyze the six filter models, and the DISC data can be 
used for two approaches to the overall system design. 

1. Select filter gains and estimate system performance given levels 
of sensor noise. 

2. Select sensor noise levels and, therefore, sensor quality and 
cost based on desired system performance. 

Once the sensor noise levels and system performance have been 
determined, the gyro filtering technique can be summarized as follows; 

1. Select filter gains from the DISC analysis based on noise and 
dynamic performance of desired '.system. 

2. Store filter models (as given in fig. 13 and figs. 17-21) and 
gains in the, VIP microprocessor memory. 

3 . Implement filter model selection algorithm in the VIP microprocessor 
based on the video sensor integration time control setting. 
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An additional consideration to the filter model selection is the 
inclusion of the drift estimator. As noted in the simulation results, 
the drift estimator causes poor disturbance torque response and slower 
response to joystick commands. It may be desirable to include the drift 
estimator only when joystick commands, or low-noise pointing, are not 
required. The drift estimate 'Could then be stored and used as a gyro 
input, but it would not be continuously updated at each computer cycle. 
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CHAPTER V. CONCLUSIONS 


Multi-star processing and gyro filtering techniques for the 
advanced VIP system have been developed and evaluated. The multi-star 
processing technique uses linearized, small-angle equations consistent 
with the limitations of the VIP microprocessor. A minimum of two guide 
stars (guide star pair) are required to uniquely detemine the three- 
axis pointing error signals. Guide star pairs and pointing error equa- 
tions must be selected based on the multi-star processing equation's 
singularities. When more than one guide star pair is available, the 
three-axis pointing error signal from each pair can be averaged to 
reduce the noise and steady-states errors in the pointing error signals. 
Consideration must be given to the number of guide star pairs used, 
since computation time increases proportionally to N, and the reduction 
in errors is proportional to the As noted in the gyro filtering 

analysis', time delays must be included in 'the filter model, and the 
noise performance (8) of the filter decreases with increased time delay. 

A further consideration for the multi-star processing is the sensitivity 
of the roll-pointing error signals as the guide stars get closer to the 
desired pointing direction. An analysis is presented that allows roll 
guide star selection based on desired roll-pointing error performance. 

The approach taken to the gyro filtering development is the use of 
a discrete, steady-state Kalman filter. A unique feature is the use of 
the integration provided by the gyro stabilization’s RIG as an integral 
part of the filter. Considerations are developed with a continuous 
filter to show how the filter can be designed and evaluated separately 
from the basic gyro stabilization. Time delays associated with the multi- 
star processing and video sensor operation are modeled in the filter as 
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unit delays of the VIP microprocessor based on the assumption of 
asynchronous operation. Six filter models and their hardware imple- 
mentations are developed. Various combinations of attitude (0), drift, 
and unit time delays are modeled. The DISC analysis can be used for 
selection of filter gains based on system performance requirements. An 
alternate use of the DISC analysis is the selection of sensor noise 
levels necessary to achieve a desired system performance; however, the 
DISC analysis strongly suggests an acceptable range of Q^/R in order 
to obtain a reasonable compromise between filter noise performance (0) 
and dynamics. 

A digital simulation of a typical gyro stabilized gimbal is developed 
and is used to validate the approach taken to the gyro filtering. Dis- 
turbance torque and joystick command responses are evaluated with two 
video sensor integration times and zero-, one-, or two-unit time delays 
included in the filters. The results confirm the use of the unit time 
delays and provide a comparative means of selecting the filter to use 
with a selected sensor integration time. A filter model with a drift 
estimate is included in the simulation, and degraded noise performance and 
dynamic response result, consistent with the DISC analysis. It is con- 
cluded that the drift estimator should not be included when precise 
joystick maneuvers or low noise pointing are required. 

Specific contributions of this work include: 

1. Development of unique multi-star processing equations which can 
be used to obtain three-axis pointing error signals from a guide star 
pair . 

2. Development'of selection criteria for equations and guide star 
pairs and for the multi— star processing that are based on equation singu- 
larities and resulting error performance. 
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3. Development of a unique, discrete Kalman filter that uses the 
stabilization gyros integration. 

4. Development of a design approach to the gyro filtering that 
separates the design aspects of the Kalman filter and the gyro 
stabilization. 

5. Development of a practical approach to the inclusion of unit 
time delays in the Kalman filter. 

While the prime objectives of this work have been accomplished, 
there are several areas where further effort could prove useful: 

1. Development and evaluation of a technique that would allow 
the measurement and use of the actual time delays in the filter 
equations at each update. 

2. Consideration of synchronous operation of the VIP microprocessor 
and the video sensor, and an analysis of the resulting performance 
tradeoffs. 

\ 

3. Consideration of a more complex method of slewing with the 
joystick during automatic pointing. For Instance, it may be possible to 
introduce slewing commands without exciting the filter dynamics by 
placing a simultaneous joystick command into the filter and at the gyro 
output into the gyro stabilization. 
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APPENDIX A 


DISC DATA FOR FILTERS 1-6 


FILTER 1 (0) 


Qd 

R 

Mq 

P0 

K 

Poles 

0.1 

0.1 

0.16 

0.06 

0.62 

0.38 

1 

0.1 

1.09 

0.09 

0.92 

0.08 

10 

0.1 

10.10 

0.10 

0.99 

0.01 

100 

0.1 

100.10 

0.10 

1.00 

O' 

0.1 

1 

0.37 

0.27 

0.27 

0.73 

1 

1 

1.62 

0.62 

0.62 

0.38 

o 

<; 

1 

10.92 

0.92 

-0.92 

0.08 

100 

1 

100.99 

0.99- 

0.99 

0.01 

0.1 

10 

1.05 

0.95 

0.09 

0.90 

1 

10 

3.70 

2.70 

0.27 

0.,73 

10 

10 

16.18 

6.18 

0.62 

0.38 

100 

10 

109.2 

9.2 

0.92 

0.08 

0.1 

100 

3.21 

3.11 

0.03 

0.97 

1 

100 

10.51 

9.51 

0.09 

0.91 

10 

100 

37.02 

27.02 

0.27 

0.73 

100 

100 

161.80 

61.80 

0.62 

0.38 
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FILTER 2 ( 0 , D ): % = 0 . 1 , T = 0.25 


Qd 

R 

<D 

Pe 

K 

Poles 

0.1 

0.1 

0.23 

0.07 

0.69 

- 0.55 

0.40 

0.77 

1 

0.1 

1.19 

0.09 

0.92 

- 0.28 

0.08 

0.93 

10 

0.1 

10.36 

0.10 

0.99 

- 0.10 

0.01 

0.98 

100 

0.1 

100.89 

0.10 

1.00 

- 0.03 

0 

0.99 

0.1 

1 

0.66 

0.40 

0.40 

- 0.25 

0.77 ± 0 . 09j 

1 

1 

1.83 

0.65 

0.65 

- 0.19 

0.38 

0.92 

10 

1 

11.22 

0.92 

0.92 

- 0.10 

0.08 
0.98 ■ 

100 

1 

101.80 

0.99 

0.99 

- 0.03 

0.01 

0.99 

0.1 

10 

2.77 

2.17 

• 0.22 
- 0.09 

0 . 88 ± 0 . 09j 

1 

10 

4.71 

3.20 

0.32 

- 0.08 

0.74 

0.92 

10 

10 

16.83 

6.27 

0.63 

- 0.06 

0.38 

0.98 

100 

10 

110.10 

9 . 17 - 

0.92 

- 0.03 

0.08 

0.99 

0.1 

100 

13.84 

12.16 

0.12 

- 0.03 

0 . 94 ± 0 . 06j 

1 

100 

17.42 

14.84 

0.15 

- 0.03 

0.92 ± 0 . 04j 

10 

100 

40.36 

28.75 

0.29 

- 0.03 

0.73 
0.98 ■ 

100 

100 

163.87 

62.10 

0.62 

- 0.02 

0.38 

0.99 
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FILTER 2 ( 0 , D ): % = 1 , T = 0.25 


Qd 

R 

M0 

P0 

K 

Poles 

0.1 

0.1 

0.38 

0.08 

0.79 

- 1.45 

0 . 42 ± 0 . 17j 

1 

0.1 

1.41 

0.09 

0.93 

- 0.81 

0.09 

0.78 

10 

0.1 

10.93 

0.10 

0.99 

- 0.30 

0.01 

0.92 

100 

0.1 

102.62 

0.10 

1.00 

- 0.10 

0 

0.98 

0,1 

1 

1.17 

0.54 

0.54 

- 0.68 

0 . 65 ± 0 . 21j 

1 

1 

2.28 

0.69 

0.69 

- 0.55 

0.40 

0.77 

10 

1 

XX • 88 

0.92 

0.92 

- 0.28 

0.08 

0.92 

100 

1 

103.56 

0.99 

0.99 

- 0.10 

0.01 

0.98 

0.1 

10 

5.07 

3.36 

0.34 

- 0.26 

0.80 ± 0 . 16j 

1 

10 

6.60 

3.98 

0.40 

- 0.25 

0-77 ± 0 . 09j 

10 

10 

18.25 

6.46 

0.65 

- 0.19 

0.38 

0.92 

100 

10 

112.16 

9.18 

0.92 

- 0,10 

0.08 

0.97 

0.1 

100 

25.34 

20.22 

0.20 

- 0.09 

0.89 ± 0 . 1j 

1 

100 

27.75 

21.72 

0.22 

- 0.09 

0 . 88 ± 0 . 09j 

10 

100 

47.10 

32.02 

0.32 

- 0.08 

0.74 

0.92 

100 

100 

168.34 

62.74 

0.63 

- 0.06 

0,38 

0.98 
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FILTER 2 (6,D): % = 0.1, T = 0.5 



R 

Me 

Pe 

K 

Poles 

0.1 

0.1 

0.30 

0.08 

0.75 

-0.50 

0.50 ± 0.05J 

V 

1 

0.1 

1.29 

0.09 

0.93 
-0.27 • 

0.08 

0.85 

10 

0.1 

10.62 

0.10 

% 

0.99 

-0.10 

0.01 

0.95 

100 

0.1 

101.69 

0.10 

1.00 

-0.03 

0 

0.98 

0.1 

1 

0.90 

0.47 

0.47 

-0.23 

0.71± 0.17J 

1 

1 

2.03 

0.67 

0.67 

-0.18 

0.39 

0.85 

10 

1 

11.52 

0.92 

0.92 

-0.09 

0.08 

0.95 

100 

1 

102.61 

0.99 

0.99 

-0.03 

0.01 

0.98 

0.1 

10 

3.93 

2.82 

0.28 

-0.09 

0.84 ± 0.13j 

1 

10 

5.62 

3.60 

0.36 

-0.08 

0.80 ± O.Olj 

10 

10 

17.49 

6.36 

0.64 

-0.06 

0.38 

0.95 

100 

10 

111.05 

9.17 

0.92 

-0.03 

0.08 

0.98 

0.1 

100 

19.80 

16.52 

.0.17 

-0.03 

0.91±0.08j 

1 

100 

22.63 

18.45 

0.19 

-0.03 

0.90±0.07j 

■ 10 

100 

43.55 

30.34 

0.30 

-0.03 

0.73 

0.95 

100 

100 

165.94 

62.40 

0.62 

-0.02 

0.38 

0.98 
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filter 2 (0,D): % = ] 


= 0.5 


0.1 0.1 0.63 0.09 

1 0.1 1.77 0.10 

10 0.1 11.80 0.10 

100 0.1 105.19 0.10 

0.1 1 1.85 0.65 

1 1 2.95 0.75 

10 1 12.88 0.93 

100 1 106.17 0.99 

0.1 10- 7.70 4.35 

1 10 9.03 4.75 

10 10 20.32 6. '70 

100 10 115.20 9.20 

0.1 100 37.40 2-7.22 

1 100 39.32 28.22 

10 100 56.18 -35.97 

100 100 174.89 63.62 


0.86 

-1.17 


0.99 

-0.30 

1.00 

- 0.10 

0.65 

-0.59 

0.75 

-0.50 

0.93 

-0.27 


0.435 

-0.24 


Poles 


0.28 ± 0.25j 


0.53 ± 0.27j 
0.50 ±,0.05j 


0.72±0.21j 

0.71±0.17j 


0.84 ± 0.13j 
0.84 ± 0.13j 
0.8 ± O.Olj 
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FILTER 3 ( 0 , 6g ): e = 0.00001 


Qd 

R 

Me 

Pe 

K 

Poles 

0.1 

0.1 

0.26 

• 0.16 

0.62 

0.62 

0 

0.38 

1 

0.1 

2.09 

1.09 

0.916 

0.916 

0 

0.08 

10 

0.1 

20.10 

10.10 

0.99 

0.95 

0 

0.01 

100 

0.1 

200.10 

100.10 

0.99 

1.44 

0 

0 

0.1 

1 

0.47 

0.37 

0.27 

0.27 

0 

0.73 

1 


2.62 

1.62 

0.62 

0.62 

0 

0.38 

10 

1 

20.91 

10.91 

0.92 

0.91 

0 

0.08 

100 

1 

200.99 

100.99 

0.99 

0.99 

0 

0.01 

0.1 

10 

1.15 

1.05 

0.1 

0.10 

0 

0.91 

1 

10 . 

4.70 

3.70 

0.27 

0.27 

0 

0.73 

10 

10 

26.18 

16.18 

0.62 

0.62 

0 

0.38 

100 

10 

209.16 

109.16 

0.92 

0.92 

0 

0.08 

0.1 

100 

3.31 

3.21 

0.03 

0.03 

0 

0.97 

1 

100 

11.51 

10.51 

0.09 

0.09 

0 

0.91 

10 

100 

47.02 

37.02 

0.27 

0.27 

0 

0.73 

100 

100 

261.80 

161.80 

0.62 

0.62 

0 

0.38 
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FILTER 4 (6,es,D) : % = J 


= 0.25, e = 0,00001 









FILTER 4 (e,eg,D); % = 1, T = 0.5, e = 0.00001 











FILTER 5 (e,6si’®s2>- ^ = 0.01 











FILTER 6 (6,0si,0s 2>D) : Djj = 1, T = 0.25, e = 0.01 


Qd 

R 

Me 

Pe 

K 

Poles 

0.1 

0.1 

2.11 

1.01 

1.51 

- 1.44 

0.42 + 0 . 17 j 





1.16 

0.01 

i 




0.79 

0.01 

1 

0.1 

6.08 

3.36 

1.34 

0.01 





- 0.81 

0.02 





1 . 15 ' 

0.06 





0.93 

0.78 

10 

0.1 

37.83 

23.^6 

1.14 

- 0.30 

0 . 00 ± O.Olj 





1.08 

0.01 





0.99 

0.92 

, 100 

0.1 

323.18 

210.28 

1.05 

- 0.00 + O.Olj 





- 0.10 

-0 





1.03 

0.98 





1.00 

' 

0.1 

1 

3.90 

2.23 

0.87 

0.64 ± 0 . 21 j 





• - 0.68 

0.01 





0.71 






0.54 

0.01 

1 

1 

7.62 

4.53 

0.97 

0.01 





- 0.55 

0.01 





0.84 

0.39 





0.70 

0.77 

10 

1 

39.07 

24.55 

1.06 

0.01 





- 0.28 

0.02 





1.00 

0.06 





0.92 

0.92 

100 

1 

324.21 

211.26 

i .04 

0 . 00 ± O.Olj 





- 0.10 

— 





1.03 

0.01 





0.99 

0.98 

0.1 

10 

10.51 

7.31 

0.46 

- 0.26 

0 . 80 ± 0 . 16 j 





0.40 

0.01 





0.34 

0.01 

1 

10 

14.44 

9.97 

0.52 

- 0.25 

0.77 + 0 . 09 j 





0.46 

0.01 





0.40 

0.01 

10 

10 

47.19 

31.73 

0.74 

0.01 





- 0.19 

0.01 





0.70 

0.37 





0.65 

0.92 

100 

10 

333.62 

220.25 

0.96 

0.01 





- 0.09 

0.02 





0.95 

0.06 




‘ 

0.92 

0.98 
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FILTER 6 (e,0gj,0g2,D): % = 1, T = 0.25, e = 0.01 

(Concluded) 
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FILTER 6 (e, 63 j,eg 2 »D): % = 1, T = 0.5, e = 0.01 










FILTER 6 (9,8sj_,0g2,D): % = 1, T = 0.5, e = 0.01 

(Concluded) 



114 








APPENDIX B 


COMPUTER SIMULATION LISTING 


Filter 4 (0,D,0g) 


a***C0NUNU0US system MODEUNG PR&^.,AM***Ar 
*****TSS/36u VERSION I MOOIFICATION 2***** 
*****t»***PROBLEM input STATi:MENTS********<» 


IMT 

THAVEeOiO 

ThCCClsytro 

- — T-I-N T -e o ^y 

thavd=o„o 

THAVDjsy^O-- 

TDCaOyO 

CriLDsa*© - - 

UCDaOyO 

TH ES T -cO: ^0 

THEGTD»yl'o 

DCSTDsO^O - - — - — 

GCMDsy^y 

TKS|C=^© . 

dynam 


DKECsINIGReIOMIC, T) 

Ta(TD»TST*TM^*57»3yjG .. - - — - 

TKsKI*ILIM 

IWM«lrWlT4It.iMM,.IUMlVD 

IaCVCPP-VB)«2*/R 

-ySaKBAUREL ^ 

URELsOMCG*»GAMO 

- TSTiC^O 

VCPPaPWAG*(VCB-VCO) 

VCOaUHmVC«#.VCP,.VCCX) . — 

VCCC«INTGRLiVCCCIC,VCC) 

ycc»i„rDi.AG(zcCyPrc,ycj 

VCsLCDLAGIZCrPCfVFOO) 

VrOOaiVFOtVfa^JkGC 

VFOauiMmvrM^vrPfVF) 

VF»CMRXPLlVF.li:ifVFlC2,LVF.lfVF2,yDi:MQG) . .. 
VCEMC0cij(|*5<i^{*25e*0G0*yF2*yF2/57?3 

O GOaOCtCFILQ 

UG»RCALPUOGlC,OCPPf OCi) 

OGXaINTCRL 40 GIlC^ts 4 - - 

GlNaOMCG4CD+GCMO+GN 
GNsGAUSSn,CYM,CSDEV) 

GDaOSO 

TDaOCli' 

JYSTaO'^o 
yCBsO^O 
VFB»C,0 
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WOG»0*0 

THE?ii»5«»^iHr*T-H£r - 

UMECA=3(jdO¥*OMEC 

THCCD»lNTCRL-f^%«, FWE-r> - - — 

PRQCEP THAV,THCCD13BE0CKA( tint, TCCD, TIME, THCCOJ 
X-F 

1 CO TO 3 

2 TNAV»(THCC-D«T«CCD4-)/TCCD- - - 

ThCCDla^HCCO 

TINT»7XWT+TCC0 - _ . . . - ... . 

CO TO 3 

. .-;_£4H4-T-X- N U E 

CNDPRO 

£CDN8CAUSS-tt-,-CCOW,-CSD£V)-- - — - 

TKAVE«ThAV+CCDN 

PROCEC THAVCTCF-i4„O,CCMD8BL0CK6t-I0C,TIMP.^ptE,F.Cl^F62,THAVi:^.ST,0G^). 
IF(TX«E-TOCF/4r5,5 ^ 

A - GQ TO -<» : t 

5 b'O.GsOG-OGO 

THAVOsTflAV^Ol - — - - 

EHRQRbJYST+THAVD-IHESTD 

- tHe-ST*5PGI-ftETW0«+^HCSTl>+X)0S-- -- 

GFIEOsTHEST-OC ■pli.TE^H 

&£ ST8ERROR * FCa>DESTD 

GCHDb-I.OaOEST 

THAVDialHAVE - - - - 

TKCSTOsTHESr 

DESTDsOgSF • - 

OGDaOG 

TPC-8TPC4.TIHP : 

GO TO 6 

6 CGWTINJJE- ..... V 

. CNDPRCr ■ 

TITLE AIROSCOPE t i -ELVATION SIMULATION. . 

TIMER DELTa«Oo5,FINriMs25,»PRDCLsiO2f0UTDEL=a(42 

■ M E THO P RK Snt 

PRINT TKETA,0MECA,GCMD,UG0,CN,THAV£,THAVD,GFIL0 

INCCN VCCCXCs««ta- - 

INCOM TlCxO^O^OMlCaO^O 

INCON VriC|50,6,VriC2sO,A . 

INCON UGICsg^OrOGIICsO'iQ 

.. p.A4UM-J£a34»» Ala> KRawOR^^ I » R a Q, 

PARAH ILIHMs-lo,,aiHP= 10 , 

PARAM VCM*-l5,,VCP=t-5T,~-- - ' - ---- - 

PARAK PWACbio< 

PARAM 2Ca-.WPCs.*.0«A74. 

PARAM ZCCa.;i,PCCabO«288 

- ..fiAR AN-CCb^^^T 

PARAK VrM*-l5,,VFpsl5, 

PARAH VFiaC39,Vr2-B37b-, ... 

PARAH 'yGHB-6*,0CPa6» 

PARAM 0GRPsti,0v25 

PARAH CCDMB0W0fCSOEVs*002788 

PARAK -PGlaJl- 

PARAM FG2B-<rt55 

P A RAH TCCD«^, 3S 

PARAH T.IMPs«fl 

END - _ 

STOP 
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